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ABSTRACT 

The recent detection of the hydrocarbon anion CeH" in the interstellar 
medium has led us to investigate the synthesis of hydrocarbon anions in a variety 
of interstellar and circumstellar environments. We find that the anion/neutral 
abundance ratio can be quite large, on the order of at least a few percent, once 
the neutral has more than five carbon atoms. Detailed modeling shows that 
the column densities of CeH" observed in IRC+ 10216 and TMC-1 can be re- 
produced. Our calculations also predict that other hydrocarbon anions, such as 
C4H~ and CgH", are viable candidates for detection in IRC+10216, TMC-1 and 
photon-dominated regions such as the Horsehead Nebula. 



Subject headings: astrochemistry — ISM: abundances — ISM: molecules — ISM: 
clouds — stars: carbon 
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Introduction 



The recent detection of CrH in the carbon-rich AGB star IRC+10216 and the cold, 



dense interstellar cloud TMC-1 by [McCarthy et al.l (120061 ) with abundance ratios relative to 
CqB. of 0.01-0.1 indicates that anions may play a more significant role in interstellar physics 
and chemistry than heretofore believed. 

The possibility that a relatively large fraction of m olecular rnaterial in interstellar clouds 
might be in the form of anions was first suggested by iHerbst I (Il98ll ) who pointed out that 
carbon chain molecules and other radicals have lar ge electron affinities, leading to high ra- 
diative attachment rates such as those measured by lWoodin et al.l (Il980l ) with the attendant 
possibility of anion/neutral fractions on the order of a few percent. More recently, it has been 
recognized that PAH anions could soak up a significant fraction of the free electrons in int er- 
stellar clouds and alter the charge balance to a significant degree (iLepp fc Dalgarn ol ll988ah as 
well as providing significan t heating thro ugh photodetachment of electrons in diffuse clouds 
( ILepp fc Dalgarno Ill988bl ). iPetrid (Il996l ) has investigated the synthesis of CN~ by the dis- 



sociat ive attachment of MgCN and MgNC and other mechanisms, while iPetrie fc Herbst 



( 119971 ) showed that C3N could be detectable in interstellar clouds. The study of la rge hy- 



drocarbons anions received some attention following the observation by iTulej et al.l ([1998) 
that absorption bands in several carbon-chain anions coincided with several of the diffuse 
int erstellar bands. Su bsequently, the formation of such species in diffuse clouds was studied 



bv lRuffle et al.l ( 



19991) wh o showed, in particular, that Cy was unlikely to be a source for any 



DIBs. iMillar et al.l (120001 ) considered the formation of hydrocarbon anions containing more 
than six carbon atoms in the carbon-rich circumstellar envelope of IRC+10216 and showed 
that appreciable column densities could arise in the outer envelope. 



Chemical Model 



The basic route to the formation of anions is electron radiative attachment: 



X + e- 



X + hu. 



which has been discussed 
The rate coefficients of 
low temperatures and in 
tachment with radiative 



in the context of bare carbon chains by iTerzieva &: HerbstI (120001 ) . 
relevant processes have not been experimentally determined at 
this work we have used phase-state theory assuming s-wave at- 
stabilisation occurring by vibrational and electronic transitions 
( ITerzieva fc Herbstll2000l ). Table [1] contains some relevant attachment rate coefficients. The 
calculations of rates for these and other processes will be discussed in a separate paper 
(Herbst, in preparation). For C„H radicals, we find that radiative attachment occurs at the 
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collisional rate once the numb er of carbon atom s is lar ger than five; for C4H, the attachment 
efficiency is only around 1%. iBarckholtz et al.l (I2OOII ) have shown experimentally that nei- 
ther carbon-chain anions nor hydrocarbon anions react with H2 but do so with H atoms. We 
have also included loss reactions of anions with C, C+, €2^2 IRC+10216), and photons, 
with all cation-anion rate coefficients taken to be 10"'' cm^ s~^. For the photodetachment 
rates, some of which are shown in Table [H we have assumed the cross section a to depend 
on photon energy e via the relation 



a 



(Jo{l- EA/ey\e> EA 



(2) 



where ctq = 1.0 x 10~^^ cm^ and EA is the electron affinity of the neutral molecule. Relevant 
electron affinities have been obtained from theoretic al and experimental results (Herbst, in 
preparation). The UV radiation field as described bvlMathis et al.l (119831) has been adopted. 
We have extended earlier mode l s such as that by iRuffie et al.l (119991 ) , which included C~ , 
n = 7-23, and by iMillar et al.l (l2000l ). which included CnH~, n = 7-23, down to anions 
containin g four carbo n atoms following the determination of the microwave spectrum of 
C4H- bv lcupta et aP J2007h . 



In addition to the reactions discussed above, we have included in the CSE and PDR 
models radiative association bet ween carbori-chain anions and neutrals leading to carbon 
chains of up to 23 carbon atoms (IMillar et al.ll2000l ). while in the dark cloud models we have 



included reactions between anions and O atoms, with a rate coefficient of 10 cm^ s 
unless specifically measured (Eichelberger et al. 2002). 



Results 



3.1. Circumstellar Envelope 



We model the car bon-ri ch AGB star I RC+10216 using the numerical model developed 
by IMillar et al.l (120001 ) and IMillar I (120031 ). They modelled the formation of anions with 
more than six carbon atoms and showed that the anion radial distribution could extend 
beyond 10^^ cm and relative to their neutral analogues, abundance ratios approaching 0.1 
(for Cf , for example) could be achieved. Some results of our present calculation in which 
we have extended our species down to C4H~ are shown in Table [2] and Figure [H Our results 
show that column densities of 1.0 x 10^^ cm~^ for C4H~ to 2.3 x 10^^ cm~^ for CiqH" are 
achievable, with anion-to-neutral column density ratios of 0.008 for C4H to around 0.3-0.4 
for CeH, CsH and CioH. Note that although the electron attachment efficiency for C4II is 
only around 1% that of CeH and CgH, C4H~ has a relatively large column density since C4H 
is more abundant. Both C4H~ and CgH^ are candidates for detection in IRC+10216. In the 



-4- 



region interior to a radial distance of 8 x 10^^ cm, the main loss of anions occurs through 
reaction with H atoms, while mutual neutralisation with C"*" dominates in (1.8 — 13) x 10^^ 
cm. Elsewhere photodetachment is the major loss route. 



McCarthy et al.l (120061 ) measured a CgH column density of 3 x 10 



cm ^, 



approxi- 
mately 60 times less than is produced in our model. Although the column density is over- 
produced, it is important to note that the column densities of the large hydrocarbon chains 
are very sensitive to the initial abundance of acetylene adopted in the model. This results 
from the fact that the most efficient method of growth is via the addition of C2 units. For 
example, if the initial abundance of C2H2 is reduced by five, the column density of CeH 
is reduced by about 45 and that of CiqH by 100. The anion column densities are not af- 
fected as much, typically reducing by about an order of magnitude, since less acetylene also 
leads to a reduction in the abundances of H atoms and in the envelope. The observed 
anion-to-neutr al ratio is 0.01-0.1 where the range is due t o the range of reported CeH col- 
umn densities (IKawaguchi et al I Il995l : iGuelin et al. 1 119971 ) , but is in reasonable agreement 
with the calculated value of 0.3, given the uncertainties in the rates of anion formation (by 
radiative electron attachment), and destruction. 

Finally we note that the anion abundances can be greater than that of free electrons 
in the region around 5 x 10^^ cm. This refiects the efficiency of carbon-chain growth and 
electron attachment. A similar result is found in dark cloud models when PAHs are included 



(ILepp fc Dalgarndll988al ) 



3.2. Dark Clouds 

Here we adopt parameters applicable to the cold dust cloud TMC-1 {n{B.2) = 2 x 10^ 
cm~^, T = 10 K, Ay = 10 mag.) and perform quasi-time-dependent modeling. The initial 
elemental abundances of C, N and O relative to H are 7.30 x 10"^, 2.14 x 10~^ and 1.76 x 10~^, 
respectively. Figure [2] shows the time-dependent evolution of the fractional abundances (with 
respect to H2) of C4H~ and CgH" and their corresponding neutrals from 10^ yr while Tabled 
and Table m present anion column densities and anion/neutral ratios respectively, at both 
early-time (3.16 x 10^ yr) and steady state (> 10*^ yr). We adopt N(H2) = 10^^ cm~^ in 
calculating the molecular column densities. 

Steady-state abundances are very low due to the incorporation of carbon into CO but 
much better agreement occurs at earlier times in the evolution. The column density of CqB.~ 
produced in our model at ear l y time (3.16 x 10^ yr), 1.35 x 10^^ cm~^, is very close to that 



observed by McCarthy et al.l (120061 ). 1 x 10 cm . We calculate the column density of 
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the neutral CeH at early time to b e 2.58 x 10_ cm~ , which is in good agreement with the 
observed value of 4.10 x 10^^ cm~^ (IBell et al. Ill999l ). The observed anion-to-neutral ratio is 
0.025, which is in very good agreement with the calculated early-time value of 0.052, given 
the uncertainties in assorted rate coefficients. 



3.3. Photon-Dominated Regions 



We have used the code developed by the Meudon group (iLe Petit et al.ll2006l ) with the 
additional hydrocarbon chemistry described above. For this model we assume a fixed tem- 
perature of 50 K and parameters representative of those in the Horsehead Nebula; i.e., G 
= 60Go, where Gq is the interstellar UV radiation field. Ay = 10 mag. to the cloud cen- 
ter, and a total hydrogen density of 2 x 10'^ cm~^. These models are steady-state, which is 
achieved rapidly due to the enhanced UV field and fast photodestruction rates. The PDR 
models show that the anionic species formed can survive exposure to a high radiation field 
present in regions such as the Horsehead Nebula because, given the high electron abun- 
dance, the formation path of electron attachment is efficient enough to overcome the rapid 
photo detechment of the fragile anions. This high electron abundance at low to intermediate 
extinction can be inferred from the C"*" abundance in Fig. [3], which also shows that significant 
abundances of the anions arise at low extinction (Ay ~ 1.5-3 mag.). At low Ay, ~ 0.5-2 
mag., the dominant loss of anions is through mutual neutralisation with C"*" (and to a lesser 
extent through reactions with atomic hydrogen). Thus, since n(C"'") = n(e) in this region, 
the anion/neutral abundance ratio depends on the ratio of the rate coefficients for electron 
attachment and mutual neutralisation. For species with more than five carbon atoms, this 
leads to abundance ratios greater than unity; for smaller species, the ratio is less than one. 
In fact, the abundances of the anions are, in some cases, larger than the neutral species. 
Relevant results are listed in Table [51 Most notably, CsH" and CiqH^ are more dominant 
than the neutral species with fractional abundances over an order of magnitude greater than 
those of CgH and CiqH. 

Since the PDR model provides abundances perpendicular to the line-of-sight for a PDR 
seen edge on, as is the case for the Horsehead Nebula, it is more appropriate to compare 
fractional abundances rather than column densities. T he peak abundance o f C4H is consis- 



tent with observations of the Horsehead region made by lTeyssier et al.l (12004 ). Although the 
relative abundance of C4H~ (3.5%) to its neutral analog is not as high as CeH" or CgH", 
the higher abundance of the neutral species allows a significantly large absolute value of this 
anion to be formed. The CeH" anion has a peak abundance of 4.7 times t hat of the neu- 



tral species. The detection of CeH in the Horsehead Nebula was reported by iTeyssier et al. 
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( 120041 ). who calculated a peak abundance relative to H2 of 1.0 x 10~^°, about one order of 
magnitude greater than calculated. It is possible that we are missing an important formation 
mechanism in this environment. One possibility is that hydrocarbon abundances are boosted 
by destruction of PAH particles, rather than being synthesized from smaller species. In any 
case, our anion/neutral ratios are independent of the mode of formation of the neutral. 



4. Conclusion 



We find that electron attachment to hydrocarbon molecules is very efficient for species 
containing more than five carbon atoms. The anions are created so efficiently that they can 
be abundant even in regions where they are destroyed rapidly by photons, such as CSEs and 
PDRs, and in dark clouds, despite the low fractional abundance of electrons. In particular, 
we find that anions such as CeH", CsH" and CiqH" can have abundances relative to their 
neutral analogs ranging from a few percent to greater than unity, while C4H~, although 
formed relatively inefficiently through electron attachment, can be observable in astronomical 
objects, given the large column densities of C4H detected. In fact, as this paper was being 
prepared, we became aware of a preprint by Cernicharo et al. in which the detection of C4ll~ 
in IRC+10216 is reported with a C4lI^/C6H^ column density ratio of 1/7, compared to our 
calculated value of 1/17. 
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Table 1. Some reactions and rate coefficients^''' 



Ri 


R2 


Pi 


P2 


a 


/9 


7 


e- 




'^4x1- 






-u.o 




e- 








o.o\j\-o ) 


-u.o 




e- 


C5H 


C5H- 




9.00{-10) 


-0.5 




e- 


C6 


06- 




1.70(-7) 


-0.5 




e- 


C6H 


C6H- 




6.00(-8) 


-0.5 




e- 


C7 


07- 




5.00(-7) 


-0.5 




e- 


C7H 


C7H- 




2.00(-7) 


-0.5 




e- 


C8 


08- 




1.70(-7) 


-0.5 




e- 


C8H 


C8H- 




6.00(-8) 


-0.5 




e- 


C9 


09- 




5.00{-7) 


-0.5 




e- 


C9H 


C9H- 




2.00(-7) 


-0.5 




e- 


CIO 


ClO- 




1.70{-7) 


-0.5 




e- 


ClOH 


ClOH- 




6.00(-8) 


-0.5 




C4H- 


PHOTON 


C4H 


e- 


1.80E-09 


0.0 


2.0 


C5- 


PHOTON 


05 


e- 


3.00E-09 


0.0 


1.5 


C5H- 


PHOTON 


C5H 


e- 


3.70E-09 


0.0 


1.5 


C6- 


PHOTON 


06 


c- 


1.30E-09 


0.0 


2.0 


C6H- 


PHOTON 


C6H 


c- 


1.50E-09 


0.0 


2.0 


C7- 


PHOTON 


07 


c- 


2.40E-09 


0.0 


1.5 


C7H- 


PHOTON 


C7H 


e- 


3.00E-09 


0.0 


1.5 


C8- 


PHOTON 


08 


e- 


l.lOE-09 


0.0 


2.0 


C8H- 


PHOTON 


C8H 


c- 


1.40E-09 


0.0 


2.0 


C9- 


PHOTON 


09 


c- 


1.60E-09 


0.0 


2.0 


C9H- 


PHOTON 


C9H 


c- 


2.00E-09 


0.0 


2.0 


ClO- 


PHOTON 


010 


e- 


l.OOE-09 


0.0 


2.0 


ClOH- 


PHOTON 


OlOH 


e- 


1.40E-09 


0.0 


2.0 



''Rate coefficients k = a(T/300)'^ exp(— 7/r) in units of cm^ 
''Herbst, in preparation. 
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Table 2. IRC+10216 hydrocarbon anion chemistry model results 



Species 


Column density 


Peak abundance 


Observed column density 




(cm~^) 


rel. to H2 


(cm^^) 


C4H 


1.3 X 10^5 


2.1 X 10^6 


2 - 9 X 10^5 (1, 2, 3, 4) 


C4H- 


1.0 X 10^3 


1.3 X 10-s 




CeH 


5.7 X 10^^ 


1.2 X 10-6 


0.3 - 3 X 10^^ (3, 4) 


CeH- 


1.7 X 10^^ 


2.6 X 10-^ 


3 X 10^2 (5) 


CgH 


2.1 X 10^^ 


4.4 X 10-^ 


5 X 10^2 (4) 


CgH- 


5.8 X 10^3 


1.4 X 10-^ 




CioH 


5.8 X 10^3 


1.3 X 10-^ 




CioH" 


2.3 X 10^3 


8.3 X 10~^ 





Not e. — References.- (1) lAyerv et al. I 
Jl993h. (ci) kawaguchi et al I Jl99,^. (A) 
McCarthy et all fl200fih . 



fll992h. (2) 



Dayal fc Sieging 



Cuelin et al. I (119971 ) . (5) 
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Table 3. TMC-1 hydrocarbon anion chemistry model results 





Early time 


Steady state 




Species 


Column density 


Column density 


Observed column density 




(cm~^) 


(cm~^) 


(cm~^) 


C4H 


4.30 X 10^3 


3.60 X 10^° 


3.4 X 10^^ (1) 


C4H- 


5.58 X 10^° 


6.99 X 10^ 




CgH 


2.58 X 10^2 


2.78 X 10^ 


4.1 X 10^2 (^2) 


CeH- 


1.35 X 10" 


2.47 X 10^ 


1.0 X 10" (3) 


CgH 


7.76 X 10" 


1.63 X 10^ 


2.2 X 10" (2) 


CgH- 


3.27 X 10^° 


6.83 X 10^ 




CioH 


8.98 X 10^° 


3.78 X 10^ 




CioH~ 


3.67 X 10^ 


2.05 X 10^ 





Note. — Refereri c es: fl ) lOuehn et al.~l (|l982h . (2) 
(3) [McCarthy et~aD J2006h . 



Bell et al 



(119991), 
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Table 4. TMC-1 hydrocarbon anion-to-neutral ratios C„H /C„H 



n Early time Steady state Observed 



4 


1.30 X 10- 


-3 


1.94 X 10' 


-3 


6 


5.23 X 10- 


-2 


8.91 X 10' 


-2 2.50 X 10-2 (1) 


8 


4.21 X 10" 


-2 


5.42 X 10' 


-2 


10 


4.10 X 10' 


-2 


5.42 X 10' 


-2 



Note. — References: (1) 



McCarthy et al, 



(120061 ). 
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Table 5. Horsehead Nebula hydrocarbon anion chemistry model results 



Species Peak abundance Obs. peak abundance 
rel. to H2 rel. to H2 



C4H 


2.4 X 10 


-9 


3.0 ± 1.2 X 10- 


' (1) 


C4H- 


8.4 X 10- 


-11 






CeH 


9.6 X 10- 


-12 


1.0 ±0.6 X 10- 


(1) 


CgH- 


4.5 X 10- 


-11 






CgH 


5.3 X 10- 


-12 






CgH- 


9.3 X 10- 


-11 






CioH 


3.8 X 10- 


-12 






CioH~ 


5.5 X 10- 


-11 







Note. — References: (1) iTeyssier et al.l (12004J ). 
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Fig. 1. — IRC+10216 circumstellar chemistry model results of hydrocarbon and anion abun- 
dances (relative to the total number density), calculated as a function of distance from the 
center of the star. 
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Fig. 2. — TMC-1 dark cloud chemistry model results of hydrocarbon and anion abundances 
(relative to H2 density), calculated as a function of time. 
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Fig. 3. — Horsehead Nebula PDR chemistry model results of hydrocarbon and anion abun- 
dances (relative to H2), calculated as a function of Ay. 



